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Abstract 

In this article, an approximate analytical solution of flow and heat transfer for a viscoelastic fluid in an axisymmetric channel 
with porous wall is presented. The solution is obtained through the use of a powerful method known as Optimal Homotopy 
Asymptotic Method (OHAM). We obtained the approximate analytical solution for dimensionless velocity and temperature 
for various parameters. The influence and effect of different parameters on dimensionless velocity, temperature, friction 
factor, and rate of heat transfer are presented graphically. We also compared our solution with those obtained by other 
methods and it is found that OHAM solution is better than the other methods considered. This shows that OHAM is reliable 
for use to solve strongly nonlinear problems in heat transfer phenomena. 
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Introduction 

In recent years a great deal of focus has been devoted to non- 
Newtonian fluid flow and heat transfer because of various 
applications in various fields of science and engineering. In this 
regard there has been increased interest in flow and heat transfer 
problem involving viscoelastic fluids. The mathematical modeling 
of channel flow of viscoelastic fluid and related heat transfer 
problems (e.g. hot rolling, extrusion of plastic) have been the focus 
of considerable research works [1-9]. Obtaining velocity and 
temperature distribution from the mathematical model in an 
efficient and reliable manner is highly important. 

The flow through a porous medium has many applications in 
science and engineering. For instance in ground water hydrology, 
reservoir engineering, petroleum engineering, chemical engineer- 
ing, chemical reactors to agriculture irrigation and drainage and 
the recovery of crude oil from the pores of the reservoir rocks [10- 
12]. 

In particular, the characteristics of heat transfer in the case of 
porous medium is of high importance due to its wide range 
applications in for example, heat exchangers, transport of heated 
or cooled fluids, micro-electronic cooling, chemical processing 
equipment, and porous burners, etc. Many researchers have 
investigated such phenomena. Example include Trimis [13] who 
has shown that porous media can be very useful within many 
applications in energy and heat-engineering and Bassam and Abu- 
Hijleh [14] who have examined heat transfer from a 2D backward 
facing step with different porous segments and analyzed the effect 
of these layers on local and overall Nusselt numbers. Nield and 
Kuznetsov [15] have analyzed the interaction of two porous layers 
with the same porosity and permeability but different thermal 



conductivity effects in forced convection and Nemoda et al. [16] 
studied a porous burner and surface burner numerically with 
different heat conductivity and power of burners. Pilevne and 
Aydin [17] have investigated forced convection in axisymmetric 
channel with different porous layers. 

Many analytical methods like homotopy perturbation method 
(HPM) [18-19], Adomian decomposition method (ADM) [20-21], 
homotopy analysis method (HAM) [22] have been successfully 
applied for different heat transfer phenomena. 

Optimal homotopy asymptotic method (OHAM) is an approx- 
imate (or semi) analytical technique that is straightforward to use 
and does not require the existence of any small or large parameter. 

The basic idea of optimal homotopy asymptotic method was 
initially introduced by Marinca and Herisanu [23]. OHAM 
reduces the size of the computational domain and can be applied 
to wide variety of problems. OHAM is a consistent analytical tool 
and it has already been successfully applied to a number of 
nonlinear differential equations arising in science and engineering. 
OHAM has been used to study steady flow of a fourth-grade fluid 
flow through a porous medium [23], oscillators with discontinuities 
and fractional-power restoring force [24], periodic solutions for the 
motion of a particle on a rotating parabola [25], thin film flow of a 
fourth-grade fluid [26], nonlinear heat transfer equations [27], and 
nonlinear problem in elasticity [28]. 

By means of OHAM, Islam et al. [29] investigated Couette and 
Poiseuille flows of a third grade fluid with heat transfer analysis, 
Idrees et al. [30] analyzed the KDV equation and Mohsen et al. 
[31] studied viscous flow in a semi porous channel with uniform 
magnetic field. Ghoreishi et al. [32] provided a comparative study 
for nth-order integro-differential equations, Mabood et al. [33,34] 
investigated heat transfer in hollow sphere and analyzed boundary 
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Figure 1. Schematic diagram of the problem. 
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layer flow, Khan et al. [35] studied thin film flow in porous 
medium. Babaelahi et al. [36] examined the viscoelastic boundary 
layer fluid flow over a stretching sheet. 

In this paper, optimal homotopy asymptotic method is utilized 
in order to study the flow and heat transfer of a viscoelastic fluid in 
an axisymmetric channel with porous wall and the influence of 
various parameters on the dimensionless velocity and temperature. 

Problem Formulation 

We consider the phenomena of flow and heat transfer of a 
viscoelastic fluid in a channel. The schematic diagram is presented 
in Fig. 1 . The x-axis is parallel to the surface of disc and z-axis is 
normal on it. The porous disc of the channel is placed at Y = +H. 
The wall that coincides with x-axis is externally heated and other 
perforated wall viscoelastic fluid is injected consistently in a way to 
keep cool the heated surface. In this regard, we follow the 
approach of [9] . 

We assume that the flow field is stagnation flow with injection. 
The flow of viscoelastic fluid is considered to be steady, 
axisymmetric and two-dimensional. The governing equations for 
the flow and heat transfer are [8]: 



(Continuity equation} 



8ru r 
dr 



dru z 

~dz~ 



Table 1. Comparison of f(rj) values using different methods. 
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The boundary conditions for Eqs. (1) - (3) are: 



u r = u z = 0, T=T W as z = 0, 
u r = 0, u z = — V, T=Tq as z = L. 



(4) 



Where u r , u z are the velocity components along in r andz — 
directions, respectively. T rr ,T (r ,T r ,-,T rz be the components of stress 
matrix, P,p,C,T,K and CD are pressure, density, specific heat, 
temperature, thermal conductivity and dissipation function 
respectively, and the dissipation function and stress components 
are defined as [8]: 



Bur u r 8u z (8u r 8u z \ 
-- Trr — + Tgg — + 1 ZZ — + z rz ( — + — )• (5) 
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% TZ = ^A rz + <j) 2 A% + <j)^B rz . 



(6) 



The similarity solution for Eqs. (1) - (3) with boundary 
conditions Eq. (4), the dimensionless similarity variables [8]: 
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Table 2. Comparison of 0(<7)values using different methods. 
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where prime denotes differentiation with respect to t]. 

Using Eq. (7) and eliminating the pressure term, Eqs. (1) - (3) 
reduced to the following ordinary differential equations or 
similarity equations. 
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where K\ = ^-,Ki = are cross viscosity parameters, Re is 
pL 1 pL 1 

the Reynolds number and Pr is the Prandtl number. Using Eq. (7) 
the boundary conditions (Eq. (4)) can be transformed for the 
considered problem are: 



/(0)=/'(0)=/(l)-l=/'(l) = 0 (10) 
0(0) = 1, 0(1)=O (11) 



The above system of Eqs. (8) - (9) with boundary conditions 
Eqs. (10) - (1 1) were studied by Kurtcebe and Erim [8] for K 2 = 0. 
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Figure 2. Effects of Aion dimensionless velocity profile for different values of Reynolds number. 
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Figure 3. Variation of dimensionless temperature with transverse distance for different values of parameters. 
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In this paper, we reconsider these equations as: 
/"" + 2Re/J'" - Ki Re(4/"/'" + 2/'/" 



9"-Pr.Re(n/'0-2/e')=O, 



(12) 
(13) 



Basic Principles of OHAM 

We review the basic principles of OHAM as expounded in [27- 
28,34—35] and other papers in the following steps: 
(i) Let us consider the following differential equation: 



A[v/(z)]+h(z)=0, tgQ, 



(14) 



and solve by means of optimal homotopy asymptotic method. 
Quantities of physical interest are the local friction 
factor, and the local Nusselt number. Physically, 

Cf r = y 2 {^Qj = ^_ _ /"(0) represents the friction factor, 
and Nu r = — — = — vRe 6'(0) defines the heat transfer rate. 



where Q is problem domain, A(y/) =Z,(w) +JV(w), where L, 
and N are linear and nonlinear operators respectively, w(t) is an 
unknown function and /z(r)is a known function, 

(ii) Construct an optimal homotopy equation as follows: 

(l-r)[L(<l>(r, r)) + h(i)] -H(r)[A(^z; r))+h(z)] = 0, (15) 




Re 



Figure 4. Variation of skin friction coefficient with Reynolds 
number and different value of K\. 
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Figure 5. Variation of heat transfer rates with Reynolds 
number and different values of Pr. 
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Figure 6. Variation of heat transfer rates with Reynolds 
number and different values of n. 
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whereO < r < 1 is an embedding parameter,i/(r) = S ]r'Ct is 

(=1 

an auxiliary function on which the convergence of the solution is 
dependent. The auxiliary function H(r) serves also to adjust the 
convergence domain as well as control the convergence region. 

(iii) If <j)(z;r,Cjj is expanded in a Taylor's series about r, the 
following an approximate solution is obtained: 



Solution of the problem via OHAM 

According to the OHAM, Eqs. (12) - (13) can be written as: 



(l-p)(f"")-H(P) 

if"" + 2Re ff" - K X Re(4/"/'" + iff"")) = 0, 



(19) 



(1 - q) (0") - G(q) (0" - Pr .Re(«/'0 - 2/0') ) = 0. (20) 



where primes denote differentiation with respect to r\. 
We consider f,9, H(p), and G(q) as follows: 



'f=fo+pfi+p 2 fi, 
= 6 0 + q6 i +q 2 6 2 , 
H(p)=pC l +p 2 C 2 , 
G(q) = qC 3 + q 2 C 4 . 



(21) 



Using Eq. (21) in Eqs. (19-20) and after simplification as well as 
rearranging the terms based on the powers oip and q, we obtain 
zeroth, first, and second order problems are: 

Zeroth order problem: 



i(2) 
0 



o ( 0 2 \>i)=o, 



(22) 



with boundary conditions: 



(v,r,Cj)=wo(x) + J2w k (z,Cj)r', j= 1,2,3,... (16) 



/o(0) = 0, / 0 (1) (0) = 0, / 0 (1) = 1, / 0 (1, (1)=0, 
0 O (O) = 1, 0 O (1)=O. 



(23) 



It has been observed by previous researchers that the 
convergence of the series Eq. (16) depends upon 
Cj , (j = 1,2,...). If it is convergent then: 



w = u' 0 (t)+ ^w*(r, Cj) 



(17) 



is obtained. 

(iv) Substituting Eq. (17) in Eq. (14), results in the following 
residual: 

R(z-Cj)=L(w(z; Cj))+h(z)+N(v/(%; Cj)). (18) 



If R(z; Cj) =0, thenw will be the exact solution although for 
nonlinear problems this will not usually be the case. For 
determining Cj , (J = 1,2,...) methods such as Galerkin's method 
or the method of least squares can be utilized. 

(v) Substitution of these constants into Eq. (17) results in the 
approximate solution. 



Its solution is: 



Mn) = i n 2 -2 n \ 

9 0 (r,) = l-r,. 



(24) 



First order problem: 

ff>(ri,Ci) = Ci (-24Re(3; 7 2 -2^) + #iRe(576 f /-288)), 
e {2) (, 7 ,C 3 ) = - C 3 Pr Re (n (6»; - 617 2 ) ( 1 - (?) - 4i/ 3 + 6; 7 2 ) , 

with boundary conditions: 



(25) 



/i(0) = 0, // 1) (0) = 0, /,(1) = 0, / 1 (1) (1)=0, 
0i(O)=O, 0i(l)=O. 



(26) 
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Its solution is: 



f {r\,C\ ) = - 24C, Re (~ ,, 7 + ~ ,« + l -Kv? - \k x ,, 4 
1 /108 288 \ , 



1 



26 



24 



(27) 
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^C 3 PrRe+ ^C 3 nPrKe ) //. 



/(Ci,C2,...,C„ ( 



where -R is the residual, 

a/ _ 3/ 
5CT ~ dCh 



Er(ri,C\,C2,...,C m ) drj, 



8J 
dC„, 



(29) 



(30) 



In case of Re = 1.5, Pr =6, K\ =0.1, n = 2, then the values of 
C\, C2, C3, C4 are given below: Ci = — 0.5997362668, C2 = 
0. 141 8766522,C 3 = 0.1 7785 14446,C 4 = 0.6537419068. 



and so on ... 

Set p=l, 9=1, to obtain the three terms solution using 
OHAM for dimensionless velocity and temperature profile as: 



f{n, C, , C 2 ) =/ 0 fa) +/, (r,, C, ) +f 2 (v,CuC 2 ), 
d(ti,C 3 ,C 4 )=9 0 (ri) + e l (r,,C 3 ) + e2(v,C 3 ,Q). 



(28) 



The method of least squares is used to obtain the convergence- 
control parameters C\, C2, C3, C4 in Eq. (28), i.e. 



Results and Discussion 

In this paper, we have successfully employed OHAM to obtain 
the approximate analytical solutions of flow and heat transfer 
problem for viscoelastic fluid in a channel with porous wall. 
Various values of different controlling parameters for both 
dimensionless velocity and temperature profiles are obtained. 
The approximate analytical solution for dimensionless velocity and 
temperature using OHAM are compared with those obtained by 
HAM [9] and numerical shooting method [9] shown in Tables 1 
and 2. These results indicate that OHAM is a feasible and efficient 
technique for problems arising in heat transfer although more 
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studies need to be conducted. Percentage of error is calculated for 
obtained OHAM solutions in Tables 1-2. The percentages reflect 
that our solution is more accurate than HAM [9] solution. 

The influence and effect of different parameters on velocity and 
temperature are displayed in Figs. 2 and 3. In Figs 2(a) and (b), the 
effect of parameters K\ and Re on the dimensionless velocity are 
displayed while keeping fixed values of parameters n and Pr. It is 
observed that the dimensionless velocity at the center of the 
channel increases with an increase raK\ . The dimensionless 
velocity also increases with Reynolds number inside the channel. 
In Fig. 3(a), the effects of Prand n on the dimensionless 
temperature are displayed for different values of K\ and Re. It 
is found that the dimensionless temperature decreases inside the 
thermal boundary layer with increasing n and Pr. The effects of K\ 
on the dimensionless temperature are found to be negligible 
whereas the dimensionless temperature decreases with an increase 
in the Reynolds number. In Fig. 4 we present and highlight the 
influence of Reynolds and cross viscosity parameter on skin- 
friction factor. It is notice that skin-friction factor coefficient 
increases linearly with an increase in Re and the opposite nature of 
skin-friction factor can be observed with the increases in cross 
viscosity parameter K\. The variation of the Nusselt number 
(representing the dimensionless heat transfer rate at the surface) is 
presented for different parameters in Figs. 5-6. Fig. 5 depicts the 
behavior of heat transfer rates against Reynolds and Prandtl 
numbers. A monotonic increase in heat transfer rate is observed 
for increasing values in both parameters. Fig. 6 shows the heat 
transfer rate against Reynolds and for different fluids according to 
power law index in). The heat transfer rate increases with the 
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